Mathematical model of vortex penetration phenomenon 
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Vortex penetration affects the stability of a superconducting system and limits the possible ap- 
plication of the system. However, the mathematical description to this phenomenon is currently 
unavailable. Here I presented a mathematical model in which I considered the effects of bulk pinning 
and internal field repulsive force on vortex hopping. Thereafter, 1 proposed a series expansion to 
the activation energy and derived a general formula to describe the time dependence of the vortex 
penetration process. With these formulas, I can analyze the experimental data and calculate the 
activation energy of the vortex penetration phenomenon. The results are accurate for the time 
dependence of the internal field measurements in a Bi2Sr2CaCu20&+ x superconductor. 

PACS numbers: 74.25.Op, 74.25.Uv, 74.25.Wx 
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I. INTRODUCTION 

An applied magnetic field can penetrate into a type-II 
superconductor in the form of quantized vortices (each 
carry a single flux quantum $ ) when the applied mag- 
netic field increases to a value above the lower critical 
field of the superconductor The vortices first pene- 
trate into the superconductor from the surfaces which 
are parallel to the applied magnetic field, and then go 
to the center of the superconductor. 7,8 Vortex motion 
is usually retarded 9 during the penetration process be- 
cause the vortices are subjected to surface attractive 
force 7 -, pinning forces, internal field repulsive force and 
damping forc o 10 ! 11 . The internal magnetic field (or pen- 
etrated magnetic field) increases gradually and finally 
reaches a saturated value. De Gennesi^ first studied the 
vortex penetration in terms of free energy. A number 
of works studied the vortex entry conditions^—, and 
other works investigated the vortex penetration into the 
bulk of superconductors over the Bean-Livingston surface 
barrier— ~—. The experimental studies of the vortex pen- 
etration were generally carried out with Hall sensors^ 
and magneto-optical imaging techniqu e 6 ^ 28 which pro- 
vides direct observations of vortex front. 

It is known that vortex entry conditions and vortex 
penetration over surface barrier were studied. How- 
ever, the Bean-Livingston surface barrier reduces expo- 
nentially with respect to the distance to the surface of 
the superconductor j 7 - The vortex motion inside a macro 
size superconductor is then mainly determined by the 
bulk pinning and internal field repulsive force. The the- 
oretical description to this phenomenon is currently un- 
available. On the other hand, the recent experimental 
observations^ 9 " have shown that the internal magnetic 
field of a high-T c superconductor is strongly time depen- 
dent after a magnetic field (below irreversible field) is 
applied to the superconductor. The time scale depends 
on the sample geometry, temperature and applied field 
magnitude. This time dependence implies that a vortex 
system is unstable during the vortex penetration process 
and a superconducting device may not work properly at 
this stage. For the purpose of application, therefore, it 



is important to construct a mathematical model which 
includes the effects of the bulk pinning and internal field 
repulsive force, and can predicts the time scale of the 
vortex penetration process. 

In this article I theoretically show that the vortex pen- 
etration process is time dependent because of the sur- 
face attractive force, bulk pinning force and internal field 
repulsive force. First, I discussed the field dependence 
of activation energy of the vortex penetration process. 
Next, I discussed the time dependence of activation en- 
ergy based on the Arrenhius relation. Finally, I derived a 
time dependent equation for the internal field in a vortex 
penetration process. 



II. MODEL 

As mentioned before, the vortices inside a type-II su- 
perconductor are subjected to various forces and vor- 
tex motion is generally retarded. The experimental 
observations^ 9 " have shown that, at lower temperatures 
and under small applied fields, the time scale of vortex 
penetration process in a superconducting crystal sample 
extends to infinity long. Therefore, it is reasonable to 
assume that the vortices are at some metastable states 
under an intermediate applied field H c \ < H a < Hi rr , 
where H c \ is lower critical field and Hi rr is irreversible 
field. In this sense, the vortex penetration process can be 
regarded as a vortices diffusion process, i.e., a process of 
vortices hopping between adjacent pinning centers. In an 
analogy to that of flux relaxation^, we can described the 
vortex penetration process with the Arrhenius relation: 
v = VQe~ Ua / kT , where vq is attempting frequency, U a is 
the activation energy of the vortices, k is the Boltzmann 
constant and T is temperature.— This relation implies 
that the activation energy U a plays an important role in 
the vortices hopping. 

In this work I am going to use some similar concepts 
and mathematical tools used in the studies of flux re- 
laxation. But it should be emphasized that the physical 
contents of these two processes are completely different. 
In vortex penetration process the vortices are pushed into 
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the superconductor by external driving force; however, in 
flux relaxation process the vortices jump out of the su- 
perconductor because of various reasons 3 -^ Let us first 
discuss the possible expression of the field dependence of 
activation energy U a in the vortex penetration process. 

A. Field dependence of activation energy 

According to Bean and Livingston^, a vortex close to 
the surface of a superconductor is subjected to attractive 
imaging force and external field repulsive force. The ad- 
dition of these two forces results in the surface barrier to 
the vortex penetration, Ubl, which is determined by the 
sample geometry and external field magnitude H a . Early 
studies have shown thal^ 

where $o is nux quantum, A is penetration depth and 
x is the distance between the vortex and sample surface. 
The first term on the right side of Eq. (JXJ) is caused by the 
attractive imaging force, and the second term on the right 
side is caused by the applied field driving force. Eq.([T]) 
shows that Ub l {x) reduces exponentially with respect to 
x, that is, x — > oo, Ubl(x) —> 0. Thus, the vortices away 
from the sample surface are almost uninfluenced by the 
surface barrier. 

On the other hand, a vortex inside a superconductor is 
subjected to bulk pinning force and internal field repul- 
sive force, which prevents the vortices motion and reduces 
the vortices hopping rate. According to the Arrhenius 
relation, this internal field repulsive force increases the 
activation energy. Therefore, the activation energy is an 
increasing function of the internal field B. Let Ub be the 
activation energy related to the bulk pinning and internal 
field repulsive force, we have dUb/dB > 0. 

The exact form of Ub(B) is usually unknown because 
of the complexity of the interaction between vortices and 
pinning centers, as well as the reaction of vortices to ex- 
ternal field driving forces. But from a pure mathematical 
consideration, we can express Ub{B) as a Taylor series of 
B (or magnetization M)M- Let us now write out the se- 
ries expression of Ub(B) explicitly, 

n 

U b (B) = U c + Y / aiB l (2) 
l=i 

where U c = U b (0), o x = U£(0), a 2 = E#(0)/2!, • ■ ■ , 

a n = ul n \o)/n\. The parameter U c is the pinning po- 
tential inside the bulk of the superconductor. The coeffi- 
cients ai (I > 2) represent the weight of the contribution 
from inelastic deformation and interaction between the 
vortices, a 2 represents the weight of the contribution from 
elastic deformation of the vortices and a\ represents the 
weight of the contribution from Lorentz forced These 
coefficients are functions of temperature T and can be 



denoted as U C (T, A, £) and ai(T, A, £). 

The total activation energy of a vortex, U a {B), is the 
summation of Ubl and Ub{B), that is 

U a (B) = U B L + U b (B) = U + J2aiB l (3) 

i=i 

where Uq = Ubl+U c is the activation energy at vanishing 
internal field, or the activation energy at time t = 0. Uq 
is determined by the sample geometry, external field and 
pinning ability. 

From Eq.(JTJ we see that the activation energy U a (B) 
is a local function, i.e., a function of position x. This is 
not explicitly shown in Eq.© because we are interested 
in the field dependence of U a in the current work. 



B. Time dependence of activation energy 

The above discussion has shown that in vortex pene- 
tration process the activation energy U a is an increasing 
function of the internal field B. It is also known that 
B increases gradually with time t. Therefore, U a is an 
increasing function of t. 

The rate of change of the internal field, dB/dt, is pro- 
portional to vortex hopping rate. On the other hand, 
B is an increasing function of t in the vortex penetra- 
tion process, i.e., dB/dt > (in relaxation dB/dt < 0). 
According to the Arrhenius relation, we have^ 3 . 

^ = Ce- u ^ kT (4) 
dt 

where C is a positive proportional constant. Eq.(j4| 
can also be derived from the diffusion equation of flux 
flow (conservation of flux).— Comparing to that of flux 
relaxation 33 , Eq.Q has differences in sign and initial 
conditions. These cause the differences between the 
mathematical models of vortex penetration and flux re- 
laxation. 

The exact solution of Eq.Q is currently unavailable, 
but we can find an approximate solution to it. Consider 
applying a magnetic field to a superconductor at time 
t = 0. The activation energy of a vortex at this moment 
is Uq, which is also the activation energy at zero internal 
field (See Eq.©). As time increases to t, the activation 
energy increases to U a . Rewrite Eq.Q and integrate it 
on both sides, 

( U \ u ^ T dU a = fc^dt (5) 

JUo JO dB 

With logarithmic accuracy 3 ^, we obtain the following 
equation: 

U a (t) = U + kTln(l+ * ) (6) 
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where r = kT / [C (dU a / dB)] is a short time scale 
parameter—. Eq.© shows that U has a strong effect 
on the time dependence of the activation energy U a (t). 



where 



w(t) = kTln 1 + 



t 



aUo/kT 



(8) 



C. Time dependence of internal field 

Our final purpose is to obtain the time dependence of 
internal field B(t). Since the vortex motion in a vor- 
tex penetration process is retarded by various forces, the 
internal field increases gradually with increasing time. 
This relation can be obtained by canceling out the ac- 
tivation energy U a from Eq.Q (field dependence of U a ) 
and Eq.© (time dependence of U a ), that is, 



w(t) 



i=i 



aiB\t) 



(7) 



To find out the time dependence of internal field B(t), 
let us now invert Eq.© by expanding B[w(t)] as a series 
of w(t), 



B[w(t)} = 



1=1 



1=1 



kTh 



t 



Mo/kT 



(9) 



Similar to earlier studies^, we have the coefficients frz's 
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where s + 2t + 3u+- ■ ■ = l — l. On considering the symme- 
try between Eq.© and Eq.©, we obtain the coefficients 



ai by doing a commutation to the coefficients bi <-> a;, 
that is 
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Eq.© describes the time dependence of internal 
field B(t). This relation can also be measured 
experimentally.™ With these relations we can calculate 
the field dependent activation energy U a (B) using the 
following procedure: first, fitting the experimental data 
with Eq.@, we obtain the fitting parameters Uq and bi 
as shown in Figure 1. Next, substituting bi into Eq. (fTT]) . 
we obtain the coefficients a; . Finally, substituting XJq and 
ai into Eq.Q, we obtain the activation energy U a (B). 

The activation energy U a (B) is a combined response to 



the internal field B. It includes the contributions from 
the surface barrier, deformation of vortices, interaction 
between vortices and possibly other unknown sources. 
Figure 1 shows that Eq.© is accurate for the experi- 
mental data of a Bi2Sr2CaCu20%j rX single crystal. 

At lower temperatures and lower applied fields, the 
inelastic deformation and interaction between vortices is 
not significant. Therefore, we can simplify the activation 
energy U a (B) in a manner analogous to that used in the 
studies of flux relaxation^. 
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FIG. 1: (Color online) Time dependence of vortex penetra- 
tion. The scattering points are the experimental data of 
a Bi2Sr2CaCu20s+x single crystal, at 25 K under an ap- 
plied magnetic field of 0.05 T. (The data is from Ref. r29T ) 
The solid black line is the theoretical fit. The fitting re- 
sults are: B(t) = biw(t) + b 2 w 2 (t) + b 3 w 3 (t), where w(t) = 
(25ft) x ln{l + t/[Te Uo ^ 25k: >]}, U = (1.27 ± 0.01) x (25ft), 
r = 0.82 ± 0.01, bi = (0.75 ± 0.01) x 10" 4 /(25fc), b 2 = 
(-0.09 ± 0.01) x 10 _4 /(25fc) 2 , b 3 = 0.00. 



First, let us consider putting ai — (I > 2) in Eq.©. 
This is equal to ignoring the inelastic deformation and 
interaction between vortices. Thus, we obtain the non- 
interacting elastic vortices which are described by the 
quadratic activation energy 



U a (B) = U + ai B + a 2 B 2 



III. DISCUSSION 

1. Measurement of U c and Ubl- — The bulk pinning 
potential U c can be obtained by measuring the time evo- 
lution of internal field B(t) at a position x » A inside 
the sample where Ubl{x) = 0. Fitting the experimental 
data B(t) with Eq.©, we have the value of U c at the po- 
sition x. In case of random pinning centers, U c is a con- 
stant and can be represented by the statistical average 
value f^N 1 / 2 , where / is the pinning force of individual 
pinning center, £ is coherence length and N is the number 
of the pinning centers in a correlation volumei 31 ' 36 . 

The surface barrier Ubl can be calculated using 
Eq. © , but we can also measure it with the help of Eq. © . 
First, using the above stated procedure to measure the 
bulk pinning potential U c . Next, fixing a sensor at a po- 
sition x' close to the sample surface and measuring B(t), 
we obtain the value of Uo(x') using Eq.©. Finally, using 
relation Ubl(x') = Uq(x') — U c , we can find out the value 
ofU BL (x'). 

If a superconductor is anisotropic, the superconduct- 
ing parameter A is anisotropic. According to Eq.©, the 
surface barrier Ubl is also anisotropic. Therefore, when 
estimating the surface barrier we must have a clear idea 
about which crystalline axis is perpendicular to the sur- 
face. 

2. Inflection point of B(t) — t curve. — As shown in 
Figure 1, the B(t) — t curve displays a concave shape 
at short time and then changes to a convex shape with 
increasing time. Let us now discuss the possible reason 
for this phenomenon. 

The fitting results indicate that it is accurate enough 
by keeping terms up to second order. Thus, we can write 
out the time dependent internal field as B(t) = biw(t) + 



(12) b 2 w 2 (t). The second derivative of B{t) is, 



Substituting Eq.© into Eq. (fT2"]) . we have (choose one 
of the solutions which is an increasing function of time) , 



1 



(13) 



where w(t) is defined by Eq.©. 



Further putting a 2 = in Eq. (fT2")l is equal to assum- 
ing that the vortices have very large elastic modulus and 
the elastic deformation can be ignored. We obtain the 
non-interacting rigid vortices which are described by the 
linear activation energy 

U a (B) = U + a x B (14) 



Substituting Eq.© into Eq. flTTl) . we have 



B{t) = —w{t) (15) 



d 2 B 
dt 2 



= 9(t) 



1 - 



2b 2 kT 



-lull 



t 



Mo/kT 



(16) 



where g(t) = 2b 2 [kT/(Te u °/ kT + 1)} 2 . Putting d 2 B/dt 2 = 
0, we obtain the inflection point time 



t* 



Te Uo/kT ]e l- 



b 1 /{2b 2 kT) 



1] 



(17) 



If t < t*, then d 2 B/dt 2 > 0; the curve is concave. If 
t > t* , then d 2 B/dt 2 < 0; the curve is convex. One can 
easily see that t* is a decreasing function of tempera- 
ture T (note that b\ and b 2 has different sign). On the 
other hand, t* is an increasing function of Uq and Uo is a 
decreasing function of the applied field H a (see Eq.©). 
Therefore, t* is a decreasing function of H a . These re- 
sults can explain the experimental observations in Ref. [2fj| 
which show that: 1. t* becomes shorter and finally goes 
to zero with increasing H a or increasing T. 2. There 
is no concave shape in the B(t) — t curves of the thin 
films, only convex shapes were observed. This is because 
the thin film samples have small surface barrier Ubl and 
consequently small Uq, which results in the unobservable 
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short t*. 

3. Effects of surface barrier on flux relaxation. — 
Eq.([T]) shows that surface barrier affects vortex penetra- 
tion. It is then natural to expect that this surface barrier 
also affects flux relaxation. 

In flux relaxation, the external field is usually zero 
and the second term in Eq.([T]) vanishes. But the first 
term in Eq.([T]) still exists. The first term is caused by 
the surface imaging force (attractive), which helps to 
draw the vortices out of the superconductor and reduces 
the activation energy of the vortices. Early studies 32 
have shown that in flux relaxation the time dependence 
of internal field is B(t) = Ya=i °i wl {t), where w(t) — 
U c — kTln (1 + t/r) and 6/ are constants. By further 
considering the surface imaging force, w(t) should be re- 
placed by 



Eq. (fT8| indicates that the surface imaging force causes 
flux relaxation theory to be a local theory. 



IV. CONCLUSION 

The general field dependent activation energy of a vor- 
tex penetration process was proposed and time depen- 
dence of internal field was derived. Although the vor- 
tex penetration is an inverse process to flux relaxation, 
there is no symmetry between the mathematical mod- 
els of these two processes because of the surface bar- 
rier. Finally, it should be mentioned that I have ignored 
the effects of damping force on the vortex penetration 
when constructing the mathematical model. The damp- 
ing force may add significant influence on the vortex pen- 
etration under large applied field where the velocity of 
vortices is large. 
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